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A 


NASA Advanced Supercomputer Facility 


* Provides massive 
computing power to all 
of NASA 

* Pleiades: 112,896 cores 

■ 7 th biggest computer in 
the world (Nov 2011) 

' Columbia: 4608 cores 

■ Formerly 2 nd biggest 
computer in the world 

* Over 1 .3 Tflops total 
compute capability 
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meline of NAS, Ames CFD, and Space Shuttle Applications 


Cray X-MP 

0.2 Gflops 

— I 

NAS Cm y 2 
Begins 2 GA°P S 


Cray Y-MP 
2.5 Gflops 


SGI Origin 2000 
128 Gflops 


SGIAltix Pleiades 
2.3 Tflops 608 Tflops 


ARC3D 



1980 


Cray C90 
15 Gflops 


SGI Origin 3800 Columbia 
1.2 Tflops 67 Tflops 
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X 


STS-1 07: Loss 

* Columbia and crew were 
lost on Feb 1 st , 2003 

• CAIB testing showed how 
a 1 .7 Ibm piece of foam 
traveling over 770 ft/sec 
could damage RCC wing 
leading edge 

■ Simulations performed at 
Ames were integral to the 
accident investigation and 
subsequent return-to- 
flight efforts 



of Columbia 
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Support for STS-107 Accident Investigation 


■ Significant improvements to fidelity of Overflow 
CFD model of SSLV 

* Steady-state simulations of many points along 
trajectory of STS-107 

“The aerodynamic loads on the bipod ramps as 
calculated by the CFD results were well within the 
design certification limits, and were a small fraction of 
the design limits at the debris-release conditions at 
MET^81.7 seconds” 

• Time-accurate 6-DOF simulations of SSLV and 
bipod-ramp foam debris using Cart3D 



Geometry Details 




New Grids with 

Rinnrf Ranvn 



Old Grid System 





Cart 3D 
Trajectory 
Results 



Outline of 
RCC panels 



Return To Flight 


* Overflow solutions of ascent 

■■ Analyze aero loads on External Tank design 
changes 

■ Provide CFD flow-fields for debris analysis 

■ Correlation of 3% Wind-Tunnel tests 

* Debris Transport Analysis 

■ Develop next generation of debris analysis 
software 

■ Develop aerodynamic models for debris 




CFD Analysis of SSLV Ascent ’ 

* Over 400 Overflow solutions run for Retum-to- 
Flight 

* New grids generated for each ascent condition 

■ 2 hours on 32 Itanium-2 CPUs 

■ 30 to 50 million grid points each 

I * Average of -1000 Itanium-2 CPU hrs / solution 

■ -20 hours of wallclock time running on 64 Itanium-2 CPUs 

IJp*' rf* Never converges to a steady-state: aft end of ET, 
attachment hardware, plumes, etc 

PP P^. Typically run for -10,000 iterations 
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Bi-pod Ramp Removal 




Alpha = -3.5 de 
Beta = -0.4 dec 
M ET = 61 sec 
Alt = 39,600 ft 


RTF Solutions 

Addition of Ice/Frost Ramps 


tach = 1.55" ^ 


ML % 
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Wind Tunnel Test Comparisons - External Tank - Phi = 203.75° 
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IA7D0B PSPBjndit'on* M.- 1 550, a = B OB'. P = B OB'. FtejmoBs 9 =2 53 xl 0°^l, 16 eleven . 1 D OB', OB eleven = 2 00" 

IA7C3A oondilnna. W_= 1.550, a = 0.G3 1 . p = fl.DQ' Reynods fl =2 50 XI D S fl!, IB eleven = 1 0.00', OB eleven = -2 00', Run - B3Q Feint = fi. LOX Roll = 15’ 
IA7D0Baondit'£:ni: M„= 1 550, a = -0 33*. £ =-0 27'. Reynold*# = 2 50x1B 6 /ft, IB eleven = 10.B0', OB eleven - -2. OB', Run =212. Point = 4, lOXHdU = B‘ 
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[A70QAFSP 

IA7QQB PSP 

1A7D0A © 
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LOo feedline bracket redesigns 


xternal Tank Redesign Assessments 


RCS Tyvek® covers 


Bipod Ramp Removal 


+Z Aerovent Modification 


Multiple ice/frost ramp redesigns 
Ascent & entry windows airloads 
Discrete airloads data book updates 
Venting database updates 
Aerothermal support 


Space Shuttle Ascent 
Debris Analysis 




Ascent Profile and 
Debris Velocities 
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Debris Transport Process Overview 
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Debris Transport 


Ballistic debris integration: 

Steady-state CFD flowfield 

Integrate motion of point-mass subject 
to drag force due to relative local wind 
vector at current location in the 
flowfield 

Neglects effect of cross-range 
dispersions due to lift 


MET = 80.0 sec 
Mach = 2.41 
Alpha = 1.68° 
Qbar = 491 psf 


mm 



Velocity, ft 1 sec 

■ 

1100 


900 


700 


500 


300 


100 


* Debris Transport software development: 

■ Developed debris-drag models using 
Cart3D 6DOF unsteady simulations 

■ Significant improvements to debris- 

trajectory computations 
Hi I mr ' 

- Wrote software for debris collision and 

proximity detection 



■ Wrote general purpose sorting and 
filtering of impact data 

■ Millions of debris trajectories have been 
computed and analyzed 
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Debris Aerodynamics Modeling 

* Debris Transport currently requires 

■ Drag model : determines impact velocity 

■ Cross-range model : determines impact locations 

- Use Cart3D CFD methods to simulate debris released in a supersonic freestream 

- Compute hundreds of 6-DOF trajectories using a Monte-Carlo approach, varying: 
Shape 

Material properties 
Initial orientation 
Initial rotation rates 

* Have developed drag and cross range models for: 

■ Foam divots 

■ Ablator material 

■ Hemisphere ice balls 

■ Bellows ice 

■ Umbilical ice 

■ Gap fillers 





Drag coefficient (S = total area) 


Foam Drag Modeling 



Drag 


Kinetic Energy 





0.75 


0.25 


1000 


100 200 300 400 500 600 700 

Non-dimensional Axial Distance (x/V 1 ') 


6- DOF Trajectories 
Static. Max. Frontal Area 
Forced-Tumble 
Static. Min Frontal Area 


Nominal Model 
Static. Max. Frontal Area 
Forced Tumble 
Static. Min. Frontal Area 
6-DOF Trajectories 


Mach number 




Foam Cross-Range Model 

• Debris can generate aerodynamic “lift” in arbitrary 
direction during trajectory (referred to as crossrange). 

• This effect is modeled in a post-processing step. 

- Crossrange cone applied to zero-lift debris trajectories 
from ballistic code to determine possible impact points. 
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Actual Divots 
Axisymmetric Frustum 

Mildly Asymmetric Frustum 


Foam Cross-Range Data 

Data from Monte- 

Carlo CFD 6-DOF ^ j= 

trajectories used to ^B j I 

develop crossrange ^a2fflE33SiaaaiP" ] T 

cone IgB [-'4 L 

Several shapes 

used to develop 11 » 

crossrange behavior IB I 

Results can be r 4 ^ 4i 

scaled to arbitrary- |H| 

sized debris ~ ■■ 

A probability can be 
assigned to any 
location within 
crossrange cone 
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Validation With Gun Development " 

Facility (GDF) Data 

* There are two aspects to the validation effort: 

■ Validate the ability of the Cart3D code to 
simulate a 6-DOF foam trajectory by direct 
comparison against range data, alidation 

■ Validate the foam drag and cross-range 
models using the range data, validation of/ 
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Side-View Cameras and Controllers 





Axial Distance (in.) 


6-DOF Method Validation 
Ames GDF ballistic data Distance vs Time 

* Mach 2.51 , 6000 g’s deceleration 




Pitch Angle (deg.) Yaw Angle (deg.) 
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Mach = 2.74 


Exp. 

6- DOF CFD 


Time (sec.) 


Shot 5, Tripped 


Shot 3, Untripped 



6-DOr . 

JVJ b rrj o cJ VMBdsBIdfl 



G Dr baJJi 
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Axial Distance (in.) 
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Drag Model Validation 

* Low oscillation trajectory - shot 2, Mach = 3.00 


Range Data 

Tumbling Frustum Model 

Nominal Frustum Model 
Static Frustum Model 



Axial Distance (in.) 


Drag Model Validation 


Medium oscillation trajectory - shot 7, Mach = 2.81 





Axial Distance (in.) 



Range Data 

Tumbling Frustum Model 

Nominal Frustum Model 
Static Frustum Model 









150 


Drag Model Validation 

* High oscillation trajectory - shot 6, Mach = 2.46 


* 




CFD (all) 

GDF Range Data (all) 
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Cross-Range Model Validation 


Cart3D 6-dof predictions bound ballistic data 

■ CFD (all) represents several hundred CFD trajectories generated 
from offset C.G. and asymmetric models 

Even mild asymmetry generates strong crossrange 
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Additional Foam Testing 


CUBRC supersonic wind- 
tunnel foam ejection tests 

DFRC FI 5 flight ejection 
tests: 38 divots 

■ 31 supersonic divots 
trimmed in high-drag 
orientation 

■ 5 subsonic divots oscillated 
or tumbled 

■ 2 divots re-contacted and 
broke apart 



f 7.4 inch divot 
Mach 3.5 
« Q 729 psf 


■ , Deceleration matches 
nominal foam drag model 



Updated Launch Commit 

Criteria 

Determine allowable ice-ball size 
on the External Tank 
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Debris-Transport Analysis Procedure 


Compute all possible ice-debris 
trajectories 

■ Release from 7600 locations 
(blue dots) 

■ 35 different masses 


Compute impact conditions 
RCC impact kinetic energy 
■ Tile damage depths 


Map a| impact data back to 1 562 
different source zones (red-grid 
cells) 


In each source zone, determine 
largest mass which does not 
exceed component capability 



DTA By the Numbers 



33 million ice-ball trajectories computed 

■ 7600 release locations 

■ 35 masses 

■ 31 flight conditions 

■ 2 ice-ball densities 

■ 2 release velocities 

10 million executions of the dprox code 

■ 1562 subset zones 

■ 35 masses 

■ 31 flight conditions 

■ 2 ice-ball densities 

■ 3 impactor targets (tile, wing LE RCC, nosecap) 

5 billion impacts evaluated 

12,000 CPU hours used 
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List of Debris Assessments 


A 



TO umbilical ice 
Orbiter/ET umbilical ice 
RCS Tyvek covers 
Tile gap-fillers 
Orbiter blankets 
Ceramic inserts 
SSME Ice 
Flow-control valves 



ET PAL ramp foam 
ET Flange foam 
ET iceballs 
ET ice/froat ramps 
ET intertank foam 
ET feed line bellows ice 
ET feed line bracket ice 

SRB Weather-seal 
SRB phenolic glass 
SRB Ablator material 
SRB viton-coated nylon 
SRB BSM RTV 




mpact Velocity, ft/sec 


mpact Angle, deg 





Example: SRB Booster Separation Motor RTV 


MET = 80.0 sec 
Mach = 2.41 
Alpha = 1.68° 
Qbar = 491 psf 
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Day of Launch Support 

* Mission Control Center, Johnson Space Center, 
Houston 

* Debris analysis team spends the hours before 
thefaunch making sure the vehicle is ready to fly 

■ Final Inspection Team 

■ Dozens of video cameras 

EJ| Looking for ice, cracks in foam, and anything unusual 

* My job includes being able to simulate potential 
debris and provide potential impact conditions 

PJTi« Execute debris analysis on NAS computers and 
produce data in less than an hour 
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Bat Debris 
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Concluding Remarks 


CFD simulations of SSLV ascent have become a 
valuable tool for the program 

Debris transport simulation has been used to 
quantify the debris environment during ascent 

Sp Helped the program focus on mitigation of the most 
dangerous debris sources 

^ Make certain that the vehicle will only launch in a safe 
configuration 



The End 



